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Within the zero-temperature hnear spin-wave theory we have investigated 
the effect of frustration and dimerization of a Heisenberg system with alter- 
nating spins Si and S2 on one- and two-dimensional lattices. The combined 
effect most visibly appears in the elementary excitation spectra. In contrast 
to the ground state energy that decreases with dimerization and increases 
with frustration, the excitation energies are shown to be suppressed in en- 
ergy by both dimerization and frustration. The threshold value of frustra- 
tion that signals a transition from a classical ferrimagnetic state to a spiral 
state, decreases with dimerization, showing that dimerization further helps 
in the phase transition. The correlation length and sublattice magnetiza- 
tion decrease with both dimerization and frustration indicating the destruc- 
tion of the long-range classical ferrimagnetic. The linear spin wave theory 
shows that in the case of a square lattice, dimerization initially opposes the 
frustration-led transition to a spiral magnetic state, but then higher mag- 
nitudes of lattice deformation facilitate the transition. It also shows that 
the transition to spiral state is inhibited in a square lattice beyond a certain 
value of dimerization. 
PACS numbers: 75.10.Jm, 75.50.Ee 
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I. Introduction 



Heisenberg ferrimagnetic spin systems, also known as alternating or mixed 
spin systems, consisting of two sublattices with spins Si and S2 of unequal 
magnitudes with a net nonzero spin per unit cell have received consider- 
able attention recently. Several theoretical studies have been carried out 
to calculate the ground-state Drooerties and the low- lying excited states of 
an alternating S1-S2 chain. Eli'ilii'iS'iilSBEl There are two low-lying ele- 
mentary excitations, the gapless ferromagnetic spin wave excitations and the 
antiferromagnetic spin waves with a gap. Dimerization of the lattice has 
been shown to lower the ground-state energy, sublattice magnetization and 
excitation energies of a ferrimagnetic system.ll2l 

Frustration due to competing antiferromagnetic second neighbor inter- 
actions can in principle destroy any LRO of the Neel type. Ivanov et a/.i 
used spin wave expansion, density matrix renormalization group, and an ex- 
act diagonalization technique to investigate the effect of weak frustration on 
the ground state energy of a quantum Heisenberg ferrimagnetic chain. They 
have identified several critical frustration parameters. The first, called ac, 
heralds a transition from the classical commensurate ferrimagnetic state to a 
spiral state. The second special point, called and termed as the disorder 
point marks the onset of incommensurate finite-range spin-spin correlations. 
The third special point, called ay, is a first order transition point from the 
long-range ordered ferrimagnetic state with total spin Sg = N{si— S2) to a 
singlet state with Sg = 0. They found that frustration causes an increase in 
the ground state energy. 

In classical ferrimagnets, one hopes to see two transition points resulting 
from frustration: one where the order is completely destroyed, and the other 
where the classical ferrimagnetic state changes into a spiral state. 

The objective of the present study is two-fold: to investigate these tran- 
sitions in a square lattice ferrimagnet, and to study the combined effect of 
dimerization and frustration. The purpose of the latter exercise would be 
to see the effect of dimerization on the transitions induced by frustration. 
For this purpose, we will study systems with (si,S2) equal to (1,^1. 
and (|, |) using a zero temperature linear spin wave (LSW) theory .Ej The 
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choice of the spins is guided by the recent assertion that the three systems 
have different predominant characters: the first has a mixed ferromagnetic 
and antiferromagnetic character, the second is more antiferromagnetic and 
the third is more ferromagnetic in character. Vahdity of the hnear spin 

wave theory will be a serious question for frustrated systems, but it has al- 
ready been argued by Ivanov et al. that the LSW theory yields satisfactory 
results, at least for small values of frustration.! 

In order to set the background for the combined effect of dimerization 
and frustration on square lattice, we shall first use the LSWT to investigate 
the effect in ferrimagnetic chains. Section II below sets up the LSW for- 
malism for a dimerized ferrimagnetic chain in the presence of frustration. It 
has already been arguedi in discussing dimerization of square lattices that 
a proper account of various possible dimerized configurations as a result of 
lattice distortions can be taken only when the nearest neighbor interaction 
is taken as J(r) ~ -. This choice made it possible to conclude that among 
the various possibilities, plaquette configuration had the lowest ground state 
energy. But this form of the nearest neighbor interaction has its own conse- 
quences, and we shall study those for the chains also in Section II in order to 
identify its particular effects to distinguish them from the effects of dimen- 
sionality of the lattice. Effects of frustration on a dimerized square lattice 
system are then studied in section III. 



II. Ferrimagnetic chains with dimerization and frus- 
tration 

Mixed spin chain systems have recently been studied extensively within 
the spin wave approximation in both undimerized0'§ili'ii'0'0lli and dimer- 
izedBu regimes. It has been conclusively shown that the linear spin wave 
theory gives excellent results for the ground state energy and magnetization 
of ferrimagnets.i We consider a chain consisting of two sublattices occupied 
by spins Si and S2 ( Si > S2) allowing for both intersublattice and intrasub- 
lattice nearest neighbor interactions Ji and J2 respectively. We choose to 
describe this system by the Hamiltonian 
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l,n • Si,„+1 + S2,n ■ S2,n+l)], (1) 

n 

where = Ji(l ± 5), and 6 is the dimerization parameter that varies from 
to 1. The total number of sites (or bonds) is 2A'^and the sum is over the 
unit cells. 



The usual boson representation of spin operators in the two sublattices is 

= {2s,-aia^y/'a^, St,^ = bi{2s,-b%y/\ 

Sin = Si- aian, = bibn - S2, (2) 

In terms of the normal mode operators 

ttfc = MfcOfc - Vkbl, (3a) 
Pk = Ukbk - Vkttl, (3b) 
the linearized Hamiltonian in Eq. ([l|) becomes 

H = eg + Yl [Ei{k)alak + E2{k)(3lf3k\ , (4) 

k 

The ground-state energy per unit cell Eg and the energies of the two excitation 
modes Ei{k) and E2{k) are given by 

e^ = C-Y.[Ai{k) + A2{k)-ak)], (5) 

k 

Ei{k) = ^{Ai{k)-A2{k)+ak)), (6) 

E^ik) = ^{A2{k)-Ai{k)+ak)). (7) 
In these equations 



6 = ^iAiik)+A2ik)f-AB^ik), (8a) 

Ai{k) = Jp S2-a si[l-cos{2k)], (8b) 

A2{k) = Jp si — a S2 [1 — cos(2A;)] , (8c) 

B{k) = v^Afc, (8d) 



Afc = Jp\/cos'^ik) + 6^sm'{k), (8e) 
C = -JpSiS2+'^{sl + sl), (8f) 
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where Jp = ^{J^ + J ) and a = ^ is the frustration parameter. 



The coefficients u{k) and v{k), constrained by the condition u'^{k) — 



V (k) = 1, are given by 
u{k) 
v{k) 



\ 



A.jk) + A^jk) + gk) 

2m 



\ 



A,{k) + A^jk) ~ i{k) 



(9a) 
(9b) 



Staggered magnetization in the two sublattices corresponding to the spins 
Si and S2, respectively, is 



Ml 
M2 



-Si- <D> 
-<D> -S2 



(10a) 
(10b) 



where < D >=< alak >=< h\hk > is the average taken in the ground state, 
which is the Neel-hke state at zero temperature: 



<D>-- 



N 



Y.v\k) 



(11) 



with k running from — | to | which is the ffist reduced Brillouin zone. 

For a two-spin system, we can think of three types of spin-spin correlation 
functions; < S^q-SI^ >, < 5'|o-5'|„ > and < 5'f q-S'|„ > . We are interested 
in the antiferromagnetic correlations which we define as 



where 



Cn — < SI q ■ 5*2 „ > < Si Q > ■ < 5*1^ > 

= -<0>2 



< O >= — ^ cos{kn) ■ u{k) ■ v{k) 



(12) 
(13) 

(14) 



and u and v are defined in Eqs.(^. Swapan et a/.B in their linear spin 
wave analysis when fit this correlation function to e~^^^ found the inverse 



correlation length ^ ^ 



Inf^. For (si,S2) 



^1, 1), this gives ^ = 1.44, 
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whereas their variational calculation gives ^ = 0.75. Others^ fit it to the 
Ornstein-Zernike form 

C{r) ~ (15) 



and found it to be 1.01. 

tcajJlfji lOLtJl^ wave theory of ferrimagnets have been discussed ear- 
none with dimerization and frustration together. 
There is a critical value of the frustration parameter a in the linear spin 
wave theory at which the energies do not remain real, signaling destruction 
of the long range order. This critical value, that we call Oc, is strongly 
5— dependent, as shown in Fig.(l). At 5 = 0, = 2(sl+s2) • -^^^ ^ (-'-'i) chain 
this is 1/3, whereas earlier DMRG results^ gave ac = 0.28. 

It is already known that with = Ji(l ± 6) the ground state energy 
decreases with S and scales as 5^.10 It however has a more interesting 
behavior with respect to a. As shown in Fig. (2), the ground state energy 
per site initially increases with a and then decreases before the long range 
order is destroyed by frustration at ac- This is true even when there is 
no dimerization, where the results agree with those of Ivanov et a/.i who 
give values only up to where the maximum occurs. The maximum shifts 
to lower values of a with S as S'^. The curves in Fig. (2) terminate at ac 
for the corresponding S. This behavior is true for all the three spin systems 
considered. 

With two atoms per unit cell, a ferrimagnet has to have two modes of 
elementary excitations. The acoustic mode is gapless {Ei{k = 0) = 0) and 
has ferromagnetic character while the optic mode E2{k) is antiferromagnetic 
and has a gap at = 0. 

Both acoustic and optic excitation mode energies decrease as 6 increases, 
as they also do when a increases. This behavior is shown in Fig. (3). The two 
excitation modes in all the three spin systems scale with 6 as 5^and linearly 
with a. There is a critical value of a at which the elementary excitation 
modes start to soften, signaling a transition from a Neel-like spin structure 
to a spiral structure.! This critical value, that we call a* and evaluate from 
the changing signs of the slopes of the dispersion curves, is different for the 
acoustic and optic modes, and in the presence of dimerization is 5-dependent: 
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'^Icoustic — 2{s(+s'^) '^p('^ ~ '^^) (16a) 

<puc = i-Jp^ - S') (16b) 

For 5 = 0, the first of these reproduces the critical value reported by 
Ivanov et aM (denoted therein as Oc). A uniform decrease of a* with 6 leads 
one to conclude that the transition to spiral spin state caused by frustration 
is facilitated by dimerization. The spin wave theory gives different behavior 
of mode softening in the two elementary excitation modes; in the case of 
the ferromagnetic mode, the mode softening starts at an a that depends 
upon the magnitudes of the two component spins while in the case of the 
antiferromagnetic mode it is uniform for all the pairs of the ferrimagnet- 
forming spins. 

The magnetization of the two sublattices, as given by Eq. (|TO|) , decreases 
in magnitude with both 6 and a as shown in Fig. (4). The decrease with a 
indicates the destruction of magnetic order. 

The spin-spin correlations decay rapidly with the spin-spin separation, 
as noted earlier also.i When fit to the Ornstein-Zernike from. Eg. ([T^) , the 
correlation length is also found to decrease with both 6 and a, as shown in 
Fig. (5). 

The case of = 

We shall now look at the same results in the event the spin-spin interaction 
amplitude is taken to vary inversely with distance between the spins. It has 
been argued earlier! that the choice of nearest neighbor interaction as J(l±5) 
did not allow taking into account the several spin-Peierls distortions possible 
in a square lattice. A more general choice we proposed was J(r) ~ ^. In 
the case of nearest neighbor coupling, this means that the amplitudes 
in Eq.(|l]) are = which approximate to the more familiar Ji(l ± 6) 
in the limit of small 6. We shall in the following describe the results for a 
ferrimagnetic chain when the proposed interaction is taken into account. 

The critical value ac now has a completely different dependence on S. 
For all the three ferrimagnetic systems, it remains constant for most of the 
5— range, but increases rapidly at higher values of 6, as shown in Fig. (6). 
This implies that in this case dimerization does not facilitate destruction of 
long range order by frustration until it is rather large in magnitude. This 
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manifests in the variation of the ground state energy with 6 and a, as shown 
in Fig. (7). Since at 5 = 0, the two cases of coincide, therefore at 5 = 
is the same for both and is 7u-^ — r- 

The difference shows up in a* also. In this case Jp = in Eq.(|l^), be- 
cause of which a* for both acoustic and optic modes becomes 5- independent. 
This contrast between the two choices of may in fact be taken to be a 
strong point that could decide, if ever tested, which of the two forms must 
be chosen to represent dimerization. 

The ground sate energy exhibits the same general 6- and a-dependence 
as in the case of the coupling J(l ± 6). Dimerization reduces the energy 
and frustration increases it, the reduction by dimerization being larger than 
the increase by frustration. The increase with frustration is again up to a 
certain value of a after which the ground state energy shows a decline before 
Oc is approached. This value of a shifts to the lower side as 6 increases, 
disappearing almost completely in the 6^1 limit as shown in Fig. (7). The 
difference is in the way in which the ground state energy scales with the two 
parameters. For a = 0, it scales with S as -jj^ , as reported earlier,! and for 
5 = it scales with a as a^'^. 

The dependence of the excitation modes on 6 and a is also quite in con- 
trast to that for the usual case of J(l ±5) coupling. Frustration continues to 
suppress the excitation mode energies. On the other hand, while for a < a*, 
an increase in 6 pushes up both the excitation energies, for a > a*, the 
acoustic and optic modes differ in their dependence on S: the acoustic mode 
is further suppressed by an increasing 6, but the optic mode is pushed up. 
This behavior is illustrated in Fig. (8). Note also that the antiferromagnetic 
mode at k = depends upon 5 but not on a. 

Sublattice magnetization and correlation lengths show the same schematic 
decreasing behavior with 6 and a as in the case of the interaction J(l ± S), 
except for the effect of the peculiar dependence of on 6 [Fig. (6)]. 

III. Frustration on a square lattice 

There are several ways in which a two-dimensional lattice can be deformed 
in dimerization. This was discussed in detail earlierB where it was shown 
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that among the various possibihties, the plaquette configuration is the lowest 
energy deformation. To arrive at this conclusion it was necessary to take the 
nearest neighbor spin-spin interactions to depend upon the spin-spin distance 
r as ^. The consequences of this for a chain have been discussed above. 
In studying the combined effect of dimerization and the competing second 
neighbor interactions on a square lattice, it becomes imperative to work with 
this form of interaction. 

Since it has already been established that among the possible deforma- 
tions of a square lattice, the one that involves two phonons, with wavevectors 
(vr, 0) and fCLvr), forming a plaquette lattice, is energetically the most favor- 
able one,iUil'0 we will restrict our investigation to this kind of deformation 
alone. 

We will write the Hamiltonian of a ferrimagnetic square lattice as a sum 
of the nearest neighbor and the next nearest neighbor (or intersublattice and 
intrasublattice nearest neighbor) parts: 



H2 
with Jx 



H = H1 + H2 

i,j A=±l 

' 51 X! "^A.A' Sl,2i,2i ■ Si 2i+A,2i+A' + ^2,21,2^ " S2,2j+A,2j+A' 

a,a'=±i 



(17) 
(18) 

(19) 



(T3A5j5 and 



J-i,-i 
1 



2(1 + 52) 



\/2(l + 6) 
1 

72(1 - 6) 



(20a) 

(20b) 
(20c) 



The linear spin wave analysis follows the same procedure as for the chain 
above. The same equations are applicable in this case, except that the various 
coefficients have now the following definitions: 
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A^{k) = 2JpS2-|{d'^(Ji,i + J-i,i) + Ci\^(Ji,i + Ji,-i)} (21a) 

A^{k) = 2JpSi-^{cf^(Ji,i + J_i,i) + Ci'i^(^i,i + ^i,-i)} (21b) 
B{k) = r(fc)Vi7^ (21c) 

C = -2JpSiS2 + ^aisl + sl), (21d) 
r(fc) = ,Jj^ {cos{k,) + cos(A;j,))' + Jl {sm{k,) + sm{ky)f (22) 

where 

Jp = (J+l + J-l)/^ = 2{l-S^) 
Jm = (J+l - J-l)/^ = 5 ■ Jp 

C^^) = 2 ■ [1 - cos{k^ + crky)] ; r = 1, 2 and cr = ±1. 

The ground state energy per site Eg defined in Eq.(^, energies of the 
two excitation modes Ei{k) in Eqs.(|) and (^, staggered magnetization Mj 
defined in Eqs.(|IO|) and correlation length defined in Eq.(|T5|) can now be 
calculated as functions of the dimerization parameter 6 and frustration pa- 
rameter a. Setting a = we reproduce the results for unfrustrated dimerized 
ferrimagnetic square lattice.!! 

The linear spin wave theory shows that, like the chain, the ground-state 
energy of a square lattice decreases with 6 and increases with a. As reported 
earlier Ji an unfrustrated ferrimagnetic square lattice has a dependence of its 

5 

ground state energy on S as ■ We now also find that ground state energy 
scales as for any fixed value of dimerization. This is true for all pairs of 
spins forming the ferrimagnet. 

The elementary excitation spectra are plotted for the system (1,|) in 
Fig. (9) along the principal symmetry directions in the irreducible Brillouin 
zone. The same schematic dispersion relations were found for the other 
two systems. The acoustic and optic modes again have ferromagnetic and 
antiferromagnetic characters respectively, and both of them are pushed up by 
dimerization and pulled down by frustration. The optic mode at k = (0, 0) 
is 6 dependent. 

As in the chains, the competing second neighbor interaction also causes 
a transition from a Neel-like state to a spiral state, indicated by softening 
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of the excitation modes, a*, the critical value at which the transition takes 
place, in the case of square lattice is also 5-dependent: 



V2(sf + si) 

This is different from the a* in chains on two counts: it is the same 
for both ferromagnetic and antiferromagnetic modes, and even when the 
interaction is J(r) = ^, it is 5-dependent. This relation also shows that for 
6 = 0, the value of a* for a square lattice is 2v^ times larger than that for 
a chain. Moreover, unlike a monotonically decreasing a* for a chain, it is 
a function that is peaked towards higher values of S, as shown in Fig. (10). 
This indicates that while the transition to spiral state in a square lattice 
is initially opposed by dimerization, it is facilitated at larger magnitudes of 
lattice deformation. This turn around in behavior occurs at 6 = 

The variation of a* and etc with 6 brings out an interesting aspect peculiar 
to a square lattice ferrimagnet. As shown in Fig. (10), there is a value of 5 
at which ac and a* are equal, and the systems which have distorted with S 
beyond this value, which we call 6*, the destruction of order occurs before 
the onset of spiral magnetic order. The value of 6* is different for different 
ferrimagnetic systems. 

The sublattice magnetization Mi decreases with both a and 5 as shown 
in Fig. (11). For a non-dimerized square lattice the magnetization has a log- 
arithmic power law scaling behavior with the frustration parameter: j^^. 
The same scaling law was found for a dimerized plaquette. 

The correlation function defined in Eq.([T5D is calculated with 

< O >= [cos{k^n^) + cos{kyny)] ■ u(k) ■ v(k) (24) 

These correlations were found to have a more rapid decay with distance 
than in a chain. The correlation length ^ in a square lattice also decreases 
with both S and a as shown in Fig. (12). There is a clear minimum in the 
correlation length at a certain a that shifts to higher values with S. 

In summary, a simple linear spin wave theory brings out quite a few new 
features in ferrimagnetic systems under the combined effects of dimerization 
and frustration. The effects in both one- and two-dimensional ferrimagnetic 
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systems are most visible in the elementary excitation spectra. Besides the 
critical value of the frustration parameter at which the long range order 
is destroyed, there is another critical value a* at which the elementary ex- 
citations undergo a mode softening, indicating a transition from a Neel-like 
to a spiral state. The LSWT shows that dimerization facilitates this transi- 
tion. Both the critical values of a are 5-dependent. While the ground state 
energy initially increases with increasing magnitude of frustration, it reaches 
a maximum and then decreases just before a reaches its critical value etc. 
Both sublattice magnetization and correlations decrease as the strength of 
dimerization and frustration increases, indicating the loss of order. The com- 
bined effects of dimerization and frustration in the case when dimerization is 
taken as ^ are quite different from the usual case Ji(l±5). The theory also 
shows that on a square lattice, dimerization initially opposes the transition 
to a spiral state, but then beyond a certain critical value 8^-, the dimerization 
parameter facilitates the transition. In the case of a square lattice ferrimag- 
net, beyond a certain value b* of the dimerization parameter, the system 
loses long range Neel-like order before going through a transition to a spiral 
state. 
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Figure captions 

Figure 1: Dependence of the critical frustration parameter etc on the 
dimerization parameter 5 for the one-dimensional spin systems (1,|), (|, 1) 
and (|, |). This is for the case when the dimerization dependence of the 
nearest neighbor interaction is taken as = Ji(l ± 5). 

Figure 2: The ground state energy Eg of the alternating spin chain (1,|) 
vs the frustration parameter a for different values of S for = Ji(l ± 5). 
The curve for each 5 terminates at the respective etc- The maximum in the 
ground state energy occurs at an a that shifts to lower values with higher 
S. In the dimer limit {6 1), the energy monotonically decreases with a, a 
feature peculiar to the combined effect of dimerization and frustration. 

Figure 3: The elementary excitation spectra for the chain (1,|) for various 
values of the frustration parameter; (a) for S — 0.0, (b) 5 — 0.4 and (c) 
6 — 0.8. This is for — Ji(l ± 6) for which the k — optic mode is 6— 
and a— independent. With nonzero a, the difference between the two modes 
is no longer constant. Softening of the two modes for a beyond a* is easily 
discernible. The other two spin systems (|, 1) and (|, |) show the same 
schematic behavior. 

Figure 4: Sublattice magnetizations for the chain (1,|) as functions of a 
and S. The curves terminate at the respective values of etc. The behavior is 
schematically the same for the other two spin systems. 

Figure 5: The variation of the correlation length ^ for the alternating 
spin chain (1,|) vs frustration parameter a, for different values of S when 

— Ji(l ± 5). The curves stop short of the respective critical value etc 
because of the strong fluctuations that ^ experiences near this point. The 
results are schematically the same for the other two spin systems. 

Figure 6: The dependence of etc on dimerization parameter when J"^ — 
for the three spin chains. This is to be seen in contrast to Fig.(l) where 
the variation is shown for the other choice of J"^. 

Figure 7: The ground state energy Eg for the alternating spin chain (1,|) 
as a function of the frustration parameter a for different values of 6 for the 
case when = j^. Like in the case of — Ji(l ± 5) the energy has 
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a maximum at an a that decreases with b. Again, in the Umit 5^1, the 
ground state energy decreases monotonically with a. The other two spin 
systems show the same schematic behavior. In accordance with the pecuhar 
dependence of etc on 5, as in Fig. (6), the curves for higher 5 can go up to 
higher values of a. 

Figure 8: The elementary excitation spectra for the chain (1,|) for vari- 
ous values of the frustration parameter and for; (a) 5 = 0.0, (b) 5 = 0.4 and 
(c) 5 = 0.8. This is for the case with — for which the optic mode 
at /c = is 5— dependent. Because ccc in this case remains constant over a 
larger range of 5, the modes in all the three figures are plotted up to just 
before the critical value a^. The mode softening is as explained in the text. 
The other spin systems have the same schematic behavior.. 

Figure 9: The elementary excitation dispersion relations of the ferrimag- 
netic system (1,|) on a square lattice. The spectra are shown for different a 
and for (a) b = 0.0, (b) b = 0.4 and (c) b = 0.8. 

Figure 10: The dependance of a* and ac on the dimerization parameter 
b for the three spin systems on square lattice. The value of b at which the 
peak occurs is independent of the spin components of a ferrimagnetic system. 
For b> b* the long range order is destroyed before the transition to a spiral 
state can take place. 

Figure 11: The a— and b— dependence of the staggered magnetization 
of a square lattice (1,|) ferrimagnet. The same schematic behavior is shown 
by the other spin systems. 

Figure 12: Correlation length ^ vs frustration parameter a for different 
values of the dimerization parameter 5 in a square lattice (1,|) ferrimagnet. 
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Fig.(3a) 
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